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Abstract 

The behavior of neutral pseudoscalar mesons 7 r° , i] and g' in hot and dense mat¬ 
ter is investigated, in the framework of the three flavor Nambu-Jona-Lasinio model. 
Three different scenarios are considered: zero density and finite temperature, zero 
temperature and finite density in a flavor asymmetric medium with and without 
strange valence quarks, and finite temperature and density. The behavior of mesons 
is analyzed in connection with possible signatures of restoration of symmetries. In 
the high density region and at zero temperature it is found that the mass of the rf 
increases, the deviation from the mass of the g being more pronounced in matter 
without strange valence quarks. 
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1 Introduction 


Different regions of the QCD phase diagram are object of interest nowadays. Un¬ 
derstanding matter under extreme conditions is relevant to understand the physics of 
heavy-ion collisions, the early universe and neutron stars. The dominant degrees of free¬ 
dom of QCD at high T and n are supposed to be partons (a quark-gluon plasma) instead 
of hadrons and restoration of symmetries is expected to occur. While the phase transi¬ 
tion at zero chemical potential and finite temperature is accepted to be second order or 
crossover, there are indications that the phase transition with finite chemical potential 
and zero temperature is first order T], Experimental and theoretical efforts have been 
done in order to explore the p — T phase boundary. Recent Lattice results indicate a 
critical ’’endpoint”, connecting the first order phase transition with the crossover region 
at T e = 160 ± 35MeV, p E = 725 ± 35MeV 0. Understanding the results of experiments 
at BNL |3j and CERN jlj provides a natural motivation for these studies. 

Much attention has been paid to the question of which symmetries are restored. In the 
limit of vanishing quark masses the QCD Lagrangian has 8 Goldstone bosons, associated 
with the dynamical breaking of chiral symmetry. In order to give a finite mass to the 
mesons, chiral symmetry is broken ab initio by giving current masses to the quarks.The 
mystery of the non-existence of the ninth Goldstone boson, predicted by the quark model, 
was solved by assuming that the QCD Lagragian has a Ua( 1) anomaly. Explicitly breaking 
the 77^(1) symmetry, for instance by instantons, has the effect of giving a mass to rf of 
about 1 GeV. So the mass of rf has a different origin than the masses of the pseudoscalar 
mesons and this meson cannot be regarded as the remanent of a Goldstone boson. 

The study of the behavior of mesons in hot and dense matter is an important issue 
since they might provide a signature of the phase transition and give indications about 
which symmetries are restored. According to Shuryak 0, there are two scenarios: only 
SU(3) chiral symmetry is restored or both, SU(3) and Ua( 1), symmetries are restored. 
The behavior of rf in medium or of related observables, like the topological susceptibility 
|6J might help to decide between these scenarios. A decrease of the rf mass in medium 
could manifest in the increase of the rf production cross section, as compared to that for 
pp collisions (7). 

Theoretical concepts on the behavior of matter at high densities have been mainly 
developed from model calculations. The Nambu-Jona-Lasinio [NJL] type models [S] have 
been extensively used over the past years to describe low energy features of hadrons and 
also to investigate restoration of chiral symmetry with temperature or density nmm 

The NJL model is an effective quark model that, besides its simplicity, has the advan¬ 
tage of incorporating important symmetries of QCD and of exhibiting also the symmetry 
breaking mechanism. Since the model has no confinement, several drawbacks are well 
known. For instance, the rf mass lies above the qq threshold. The model describes this 
meson as qq resonance, which would have the unphysical decay in qq pairs, and the def¬ 
inition of its mass is unsatisfactory. For this reason, it has been proposed that in order 
to investigate the possible restoration of Ua{ 1) symmetry, instead of m v >, the topological 
susceptibility y, a more reliable quantity in this model, should be used [Hj. 
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The behavior of 577(3) pseudoscalar mesons, in NJL models, with temperature has 
been studied in HmHnzi' Different studies have been devoted to the behavior of pions 
and kaons at finite density in flavor symmetric or asymmetric matter mrnmrn- 
This paper is devoted to investigate the phase transition in hot and dense flavor 
asymmetric matter and the in medium behavior of the neutral 577(3) pseudoscalar mesons, 
in the framework of the SU (3) NJL model including the’t Hooft interaction, which breaks 
the 17a(1) symmetry. We choose these mesons because they might carry information about 
the restoration of symmetries: pions, which have been studied extensively in this concern, 
are a privileged tool to appreciate the restoration of chiral symmetry, since, being the 
lightest pseudoscalar mesons, are the best Goldstone-like bosons. On the other side, the 
in-medium behavior of // and r\ is expected to carry information about the restoration of 
£ 74 ( 1 ) symmetry. At variance with the other mesons, there is not much information about 
the behavior of these mesons in medium. Concerning the r/ meson, the recent discovery 
of mesic atoms might provide useful information in this concern [151. 


2 Model and Formalism 


We use a version of the 517(3) NJL model described by the Lagrangian: 


L = 


q (id • 7 - m) q + y J2 [(qX a qf + {q(i^)X a qf 


a=0 


9d 


det[?(l - 75 ) 9 ] + det[g(l - y 5 )g] 


( 1 ) 


Here q = (u, d, s ) is the quark field with three flavors, Nf = 3, and three colors, N c = 3. A a 
are the Gell - Mann matrices, a = 0,1,..., 8, A 0 = yj 11. The explicit symmetry breaking 
part ([IJ contains the current quark masses rh = diag(m u , m^, m s ). The last term in dTJ) is 
the lowest six-quark dimensional operator and it has the 517^(3) ® 5 E/r( 3) invariance but 
breaks the £ 74 ( 1 ) symmetry. This term is a reflection of the axial anomaly in QCD. For 
the general reviews on the three flavor version of the NJL model see Refs, in in la uni 
To obtain a four quark interaction from the six quark interaction we make a shift 
( qX a q ) —> ( qX a q)+ < qX a q >, where < qX a q > is the vacuum expectation value, and 
contract one bilinear ( qX a q ). The hadronization procedure can be done by the integration 
over quark fields in the functional integral, leading to the following effective action: 


w.„[vM = - 1 - - \ 


2 

-iTi In 


i{lM - rh + a a X a + {i^){^ a X° 


( 2 ) 


The fields o a and ip a are scalar and pseudoscalar meson nonets. We have introduced 
projectors 


Sab dS^ab T QD^abc ^ qX q , (3) 

Pab 9S$ab QD^abc ^ qX q 7* . (4) 
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The constants D abc are such that they coincide with the SU( 3) structure constants 
d abc for a, 6, c = (1, 2,..., 8) and D 0ab = ~^ ab , D 000 = ^/|. 

A straightforward generalization of the model for finite temperature and density can 
be done by using the Matsubara technique (see Ham). The first variation of the action 
0 leads to the well known gap equations. By expanding the effective action © over 
meson fields and keeping the pseudoscalar mesons only, we have the effective action: 



— (p u 

2 


p * ~ IMP) 




1 

2 


<p a D£{P)<p b 


(5) 


where D al f(P) is the inverse unnormalized meson propagator and fI a& (P) is the polariza¬ 
tion operator, which in the momentum space has the form 


bU(P) = iN c 


d A p 

( 2 W 


tr D 


SiW(A a )«(i7 5 )Si(p + ^)(A%(i7 5 ) 


( 6 ) 


The model is fixed by coupling constants gs,9D, the cutoff parameter A which regu¬ 
larizes momentum space integrals and current quark masses. We use the parameter set 
rn u = m d = 5.5 MeV, m s = 140.7 MeV, g s A 2 = 3.67, g D A 5 = -12.36 and A = 602.3 MeV, 
that has been determined by fixing the conditions: M n o = 135.0 MeV, Mk = 497.7 MeV, 
/,r = 92.4 MeV and Mp = 960.8 MeV. We also have: M v = 514.8 MeV, 9(Mff) = —5.8°, 
g V uu = 2.29, g V ss = -3.71 = 960.8 MeV, 9(M'*,) = -43.6°, g r fu U = 13.4, g rrss = -6.72. 

Although, due to the lack of confinement in the model, the rf - meson is here described 
as a resonance state of qq, we use as an input parameter. 

To consider the diagonal mesons 7T°, rj and g' we take into account the matrix structure 
of the propagator in (JSJ) (See Ref. jH, TTHTTIj for details). In the basis of 7T° — rj — g' system 
we write the projector P ab and the polarization operator P ab as matrices: 



( P 3 3 

^30 

P 3 8 \ 

1 n 33 

n 30 

n 38 \ 

Pab = 

Pq3 

Poo 

P 08 and U ab = 

n 03 

noo 

n 08 


\ -^83 

P80 

Pss ) 

v n 83 

n 80 

n 88 / 


It should be noticed that in media with different densities of u and d quarks the non¬ 
diagonal matrix elements P 30 , P 38 oc (<< q u q u » — « qsu >>) and n 30 , n 38 oc 
(j«u(p) _ Pjd^pyj are non vanishing and correspond to 7T° — g and 7T° — rf mixing (here 
<< q^i » are the in medium quark condensates and Pf‘(P) are write below). Here we 
did the approximation of neglecting these elements, whose effects we found to be negligi¬ 
ble by means of a simple estimation. The nonvanisnhing elements of n a & depend on the 
regularized integrals [TH] : 


J“(P 0 ,r,JUi) = 



( 8 ) 


where nf = are the Fermi distribution functions of the negative (positive) 

energy state of the ith quark. Notice that, as the temperature and/or density increase 
these integrals may acquire an imaginary part, as long as the meson masses cross the qq 
threshold. This is also true for rf in the vacuum, as already mentioned. 
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Figure 1: Pression as function of density and temperature. 

3 Discussion and conclusions 

We present our results for the nature of the phase transition and for the masses of 
7T° , p and p' in three scenarios: zero density and finite temperature, zero temperature and 
finite density and finite temperature and density. Two types of asymmetric quark matter 
that are considered: for the first case (Case I), that might be formed temporarily in 
heavy-ion collisions, we fix densities by pd = 2p u and p s = 0; for the second one (Case II), 
that might exist in neutron stars, we impose the condition of (3 equilibrium and charge 
neutrality through the following constraints, respectively on the chemical potentials and 
densities of quarks and electrons Pd = Ps = Pu + p e and | p u — \ (pd + Ps) — p e = 0, with 

P. = MP - M I? l2 0(l4 - Mf ) and p e = |h 

The nature of the phase transition in NJL models at finite T and/or p has been 
discussed by different authors nnunumnnuisi - At zero density the phase transition is a 
smooth crossover one; at non zero densities different situations occur. Here we analyze this 
question by examining the curve of the pressure as a function of density and temperature 
(Fig. 1), defined as P(p,T) = — [il(p, T) — f](0,T)], where fi(p, T) is the thermodynamic 
potential muni i^di. At zero temperature and finite density, as it has been shown in 
mm, there is a region of densities where the pressure is negative, meaning that a phase 
of low density and broken symmetry coexists with a phase of high density and partially 
restored chiral symmetry (see Fig. 1). A suitable choice of parameters leads to the first 
zero of the pressure at p ~ 0 and the second zero, (at p c = 2.25, for Case I and Case 
II ) corresponds to a minimum of the energy per particle, which means that a stable a 
hadronic phase exists at these densities, consisting of droplets of quarks surrounded by 
a non trivial vacuum. As pointed out in m different features occur in the quark phase 
( for p > p c ) according to we consider or not matter in beta equilibrium. In the last 
case the mass of the strange quark becomes lower than the chemical potential at densities 
above ~ 3.8po what implies the occurrence of strange quarks in this regime. As it will be 
discussed below, this fact leads to meaningful differences in the behavior of p , rj. 
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Figure 2: p, p' and 7r° masses as function of temperature a) and density: b) Case I; c) 
Case II. 

Now we discuss the nature of the phase transition at finite density and temperature by 
studying the pressure in the T—p that is shown in Fig. 1. One can see that there is a region 
of the surface with negative curvature that corresponds to the region of temperatures and 
densities where the phase transition is first order, since the pressure and/or compressibility 
are negative. The point T = 56Mev, p = 1.53p, were the pressure is already positive and 
the compressibility has only one zero, is identified, as usual, as the critical endpoint, that 
connects the first order and the second order phase transition regions. Above that point, 
we have a smooth phase transition, corresponding to the regions of positive curvature of 
the surface. 

Unlike at T — 0, we can not talk anymore about a stable hadronic phase in the first 
order phase transition region because the absolute minimum of the energy per particle is 
now at zero density. In spite of these drawbacks, we think it is illustrative to plot the 
masses of 7T° , p and p' in the whole region of the p — T plane. 

Let us now analyze the in medium behavior of the mesons under study. To understand 
our results, it is useful to study the limits of the Dirac sea (denoted by Ui) of this mesons. 
They can be obtained by looking the limits of the regions of poles in the integrals I%(Pq) 
(see 0): tu u = wd = 2 p d and u s = 2p s . (At finite temperature, we will have 
uj u = 2 M u , uj d = 2M d and u s = 2M s ). 

The behavior of these mesons at p = 0 and I/O has been study by (5} [T?]. We 
reproduce these results for the sake of comparison with our new findings at T = 0, p ^ 0 
(see Fig. 2a)). The mesonic masses change smoothly and when they exceed the sum of 
the masses of the constituent quarks the mesons became unbound (Mott temperature). 
Usually, the critical temperature is defined as the Mott transition temperature for the 
pion. The masses of p and 7 / show a tendency to become degenerated, even after the 
critical temperature. The decrease of the mass of p' is not generally considered enough 
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Figure 3: Mixing angle, 9, as function of density. 


to give an indication of restoration of Ua( 1) symmetry in hot media [6, [7j. 

Now, we analyze the results for the masses of p and rj at finite densities and zero 
temperature, that are plotted in Fig. 2 b)-c). Up to the critical density, they exhibit 
a tendency to became degenerated, but, after that point they split again, the splitting 
being more pronounced in Case I. These results show that if a degeneracy of 77,77' is an 
indication of the restoration of Ua{ 1) this is unlikely at high densities. In Case II the 
system shows a tendency to the restoration of flavor symmetry, which is related with 
the presence of strange quarks in the medium that occurs at about p ~ 3.8p 0 and, as 
shown in Hi, manifests itself in other observables. As we can see, the rf - meson lies 
above the quark - antiquark threshold for p n < 2.5po and it is a resonant state. After 
that density, the 7 / becomes a bound state. We checked that the behavior in matter with 
Pu = Pd, Ps = 0 is qualitatively similar to Case I. 

In order to understand the difference between the results in the high density region 
for Cases I and II, it is convenient to study the in medium behavior of the mixing angle, 
9, that is plotted in Fig. 3. As it is well known, the quark content of p and p' depend on 
the mixing angle in the form: 


I P > 
I rf > 


cos 9 


v/6 

1 


\uu + dd — 2 ss > — sin 6—=\uu + dd + ss >, 

Vo 


cos#— i=\uu + dd + ss > + smQ—=\uu + dd — 2ss > . 

V6 


( 9 ) 

( 10 ) 


Above p ~ 3.5po the angle becomes positive and increases rapidly and the strange 
quark content, y = {uu + dd)/ss, of the mesons changes: at low density, the p' is more 
strange than the p but the opposite occurs at high density. Since in Case I there are 
no strange quarks in the medium (what implies that the strange quark mass is almost 
unaffected), the p mass should stay constant in the region of densities where its content 
is dominated by the strange quark. This explains the larger splitting between p and p' in 
Case I. 


7 









Figure 4: Mesons masses as functions of temperature and density (Case II). Upper panel: 
r) and rf mesons. Lower panel: n° meson. 


Considering now the finite density and temperature case (we plot the masses of the 
mesons for the case of quark matter with strange quarks in (3 equilibrium at finite T in 
Fig. 4), and similarly to the finite temperature and zero density case, the qq threshold for 
the different mesons is again at the sum of the constituent quark masses, so the mesons 
dissociate at densities and temperatures close to the critical ones. A second feature to be 
noticed is that the mesons that have a Goldstone boson like nature show more clearly the 
difference between the chiral symmetric and asymmetric phase. The diagram for 7r° shows 
clearly a ” line” that separates the chiral broken phase (region of negative curvature) from 
the chiral restored phase. This is not so evident for the rj. A more smooth behavior is 
found for rj' (Fig. 4), which is natural, since it is not a Goldstone boson associated with 
chiral symmetry. 
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